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This  technical  note  describes  the  low-noise,  high-resolution  subbottom 
imaging  system  developed  to  remotely  and  efficiently  determine  character¬ 
istics  of  subbottom  marine  sediments  as  they  relate  to  dredging.  The  theo¬ 
retical  foundation  of  this  approach,  which  is  based  on  acoustic  impedance, 
is  described.  A  case  history  survey  of  Galveston  Ship  Channel  is  also 
summarized. 

Background 

The  focus  of  this  work,  conducted  by  the  Dredging  Research  Program,  is 
toward  developing  a  technique  to  remotely  and  efficiently  determine  char¬ 
acteristics  of  subbottom  marine  sediments  as  they  relate  to  dredging.  A 
low-noisc,  high-resolution  subbottom  imaging  system  is  essential  to  this 
program.  To  fulfill  this  requirement,  a  digital  data  acquisition  system  has 
been  combined  with  specialized  processing  software  to  accurately  assess 
bottom  and  subbottom  in  situ  conditions. 

These  undertakings  respond  to  the  fact  that,  each  year,  the  U.S.  Army 
Corps  of  Engineers  spends  millions  of  dollars  on  river  and  harbor  mainte¬ 
nance  and  ship  channel  realignment  projects.  Currently,  the  Corps  relies 
on  drilling  and  laboratory  testing  programs  to  assess  marine  sediments  in 
terms  of  material  type,  density,  and  thickness  for  purposes  of  characteriz¬ 
ing  proposed  dredging  sites.  However,  sampling  and  coring  programs  are 
costly,  provide  only  discontinuous  information  about  material  characteris¬ 
tics,  and  cannot  effectively  address  situations  where  actual  subbottom 
conditions  are  highly  variable. 

Additional  Information 

This  technical  note  was  written  by  Mr.  R.  F.  Ballard,  Jr.,  Mr.  K.  J. 
Sjostrom,  Mr.  R.  G.  McGee,  and  Mr.  R.  L.  Leist.  For  additional 
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information,  contact  Mr.  Ballard,  (601)  634-2201,  or  the  manager  of  the 
Dredging  Research  Program,  Mr.  E.  Clark  McNair,  Jr.,  (601)  634-2070. 

Theoretical  Foundation 

The  acoustic  impedance  method  is  a  modification  of  the  seismic  reflec¬ 
tion  technique  commonly  used  in  offshore  oil  exploration  but  tailored  to 
shallow-water  environments.  As  energy  generated  from  an  acoustic 
source  (in  the  form  of  a  plane  wave)  arrives  at  a  boundary  between  two 
layers  of  differing  material  properties,  part  of  the  energy  will  be  reflected 
back  toward  the  surface  and  part  transmitted  downward.  Portions  of  the 
transmitted  energy  will  undergo  absorption  or  attenuation  in  the  layer 
while  the  remainder  propagates  through  to  the  next  stratigraphic  bound¬ 
ary.  Ratios  between  transmitted  and  reflected  energy,  called  reflection  co¬ 
efficients,  are  dependent  on  the  density  and  velocity  of  the  materials 
through  which  the  energy  is  propagating. 

Wave  velocities  are  controlled  by  clastic  properties  of  the  two-phase  sed¬ 
iment  mass  (sea  water  in  pores  and  mineral  structure).  Properties  such  as 
porosity  and  grain  size  affect  sound  velocity  only  through  their  effects  on 
the  elasticity  of  the  sediment.  In  previous  studies  (Hamilton  1970,  1972), 
it  was  concluded  that  elastic  properties  of  water-saturated  sediment  could 
be  expressed  through  Hookean  elastic  equations,  unless  attenuation  is  con¬ 
sidered,  in  which  case  linear  viscoelastic  equations  are  recommended. 

The  basic  equation  for  the  velocity  of  a  compressional  wave  Vp  is 

V„  -  [(*  <  4/nMl/pl*  (1) 

where 

k  =  incompressibility  or  bulk  modulus  and  equals  (1 /p) 

P  =  compressibility 
p  =  shear  (rigidity)  modulus 
p  =  saturated  bulk  density 

When  a  medium  lacks  rigidity.  Equation  1  becomes 

V9  =  (K/p)14  (2a) 

or  , 

=  (1/pp)'4  (2b) 
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Compressibility  p  and  density  p  in  Equation  2b  have  been  expanded  into 
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where  Tl  is  the  volume  of  pore  space  occupied  by  water  (fractional  porosity) 
and  subscripts  s  and  xv  indicate  mineral  solids  and  water. 

The  influencing  parameters  of  this  basic  seismic  wave  equation  suffice 
to  answer  the  question  Why  acoustics  to  characterize  hGttomjsubboiiom  materi¬ 
als?  To  continue,  the  acoustic  reflection  coefficient  (R)  is  defined  as 


where 


Er  =  reflected  energy 

£/=  total  energy  incident  to  the  boundary  (see  Figure  1) 
The  reflection  coefficient  is  also  equal  to 

p  =  %lizA 

(Z  s  +  Z  ) 


where 


Z-w  =  P w  Cw  =  water  impedance 
Zs  -  ps  Cs  =  soil  impedance 
Pw  =  1  g./cm3 

C_  1  cn  /Vtn  — /  _  _  , 

to  —  i3u,uuu  cui/seC 


Hence,  it  is  clear  that  the  acoustic  impedance  (Zs)  of  the  surficial  layer 
can  oe  calculated  readily.  The  product  of  transmission  velocity  and  den¬ 
sity  of  material  is  the  acoustic  impedance  and  represents  the  influence  of 
the  material's  characteristics  on  reflected  and  transmitted  wave  energy. 
The  relahonship  between  acoustic  impedance  and  specific  soil  properties 
has  been  empirically  based  on  an  extensive  database  of  world  averages  of 
impedance  versus  sediment  characteristics  (Hamilton  1970, 1972;  Hamilton 
and  Bachman  1982)  (see  Table  1). 
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Table  1.  Soil  Classification  Versus  Acoustic  Impedance  Range 

(Hamilton  1972) 


Description 

Water 
Silty  clay 
Clayey  silt 
Silty  sand 
Very  fine  sand 
Fine  sand 
Medium  sand 
Coarse  sand 
Gravelly  sand 
Sandy  gravel 

Note:  Values  corrected  for  temperature  and  salinity. 


At  this  point,  it  must  be  emphasized  that  Hamilton's  pioneering  efforts 
were  limited  to  surficial  bottom  materials.  To  extend  the  depth  of  inves¬ 
tigation  into  multilayer  subbottom  environments,  Caulfield  and  Yim  (1983) 
devised  a  model  correlating  to  Hamilton's  work.  The  model  is  used  to 
correct  for  absorption  and  other  losses  in  bottom  sediments  as  a  function 
of  frequency  so  that  the  reflection  coefficients  and  acoustic  impedancp  of 
sediments  can  be  calculated  as  if  they  were  surficial  sediments  (reflectors). 
In  practice,  the  concept  is  extended  to  each  subsequent  layer  until  the 
signal-to-noise  ratio  is  at  a  level  from  which  information  cannot  be  ex¬ 
tracted  with  accuracy  (normally  5  db).  The  model  is  then  combined  with 
classical  multilayer  reflection  mathematics  to  yield  reflection  coefficients 
equivalent  to  surficial  sediments  for  subbottom  layers.  Since  some  as¬ 
sumptions  must  be  made  regarding  attenuation  factors  (determined  from 
site-specific  borings  and  laboratory  data),  this  approach  must,  be  defined  as 
empirical. 

Equipment 

A  seismic  source  of  known  energy  content  as  a  function  of  frequency, 
deployed  just  below  the  water  surface,  generates  acoustic  waves  that  prop¬ 
agate  downward  through  the  water  column  and  sediments.  High- 
resolution  profiling  Systems  specifically  designed  for  shallow-water  use 
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Acoustic  Impedance,  — r  — 

cmvsec 

1,450 

2,016-2,460 

2,460-2,864 

2,864-3,052 

3,052-3,219 

3,219-3,281 

3,281-3,492 

3,492-3,647 

3,647-3,880 

3,880-3,927 


and  operating  at  frequencies  below  12  kHz  are  typically  used.  As  a  rule, 
lower  operating  frequencies  allow  greater  energy  penetration  into  the  sub- 
bottom,  but  because  of  longer  wavelengths,  lack  the  vertical  resolution  of 
higher  frequency  systems. 

Two  commercially  available  geophysical  instruments,  a  3.5-kHz 
"pinger"  system  and  an  integrated,  high-definition  400-Hz  to  5.0-kHz 
"boomer"  system,  are  normally  used  to  fulfill  the  above  criteria.  Recently, 
tests  with  chirpers  have  been  conducted  to  assess  potential  advantages 
However,  data  discussed  herein  were  obtained  with  the  more  coriven- 
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seismic  sources  propagates  through  sediment  of  varying  densities  and 
acoustic  velocities,  energy  is  reflected  at  geologic  boundaries  where  there 
is  a  distinct  contrast  in  the  acoustic  impedance  between  layers.  Reflected 
signals  are  amplified,  filtered,  and  recorded  with  a  specially  designed  shal¬ 
low  seismic,  digital  data  acquisition  system  developed  in  conjunction  with 
Caulfield  Engineering  (Caulfield  1991a).  Energy  loss  as  a  function  of  fre¬ 
quency  is  then  determined.  Tire  system  also  provides  real-time  presenta¬ 
tion  of  the  seismic  signal  for  acquisition  quality  control. 

Data  Acquired  and  Ground  Truth  Comparisons 

Because  of  the  nonuniqueness  of  seismic  reflection  signatures,  several 
combinations  of  geologic  conditions  could  conceivably  yield  similar  signal 

:  O.l'l  -1  T*0  /it/trin  t-l  «'-*/->  M/vwt  v«(  •  (-/->  s-J  awtw  .r*  J  ~  1  ..  -  -i.  —  _ _ 1  _  .!  . 
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regions  such  as  the  Mississippi  Sound,  Savannah  Ship  Channel,  or  San 
Francisco  Bay,  differing  sediment  units  usually  have  a  characteristic  and 
relatively  narrow  range  of  impedance  values.  Therefore,  using  calibration 
procedures  that  incorporate  local  core  and  laboratory  data,  seismic  reflec¬ 
tion  data  are  processed  at  known  sample  locations  to  yield  acoustic  imped¬ 
ance  values  of  the  known  reflection  horizons. 

Estimates  of  in  situ  density  are  derived  from  computed  impedance  val¬ 
ues  and  correlated  with  ground  truth  information.  Acoustic  predictions 
versus  core  data  for  consolidated  materials  in  Mobile  and  Gulfport  Ship 
Channels  was  presented  by  Ballard,  McGee,  and  Whalin  (1992).  Results 
documented  were  within  1  percent. 

Similar  findings  have  also  been  observed  when  comparing  acoustic  im¬ 
pedance  density  predictions  to  nuclear  densitometer  data  for  fluff/fluid 
mud  type  materials  in  the  Gulfport  Ship  Channel.  Again,  correlation  has 
been  excellent.  A  continuing  program  of  database  expansion,  coupled 
with  ground  truth  information  obtained  using  a  wide  variety  of  condi¬ 
tions,  is  necessary.  Testing  to  date  has  shown  that  density  estimates  to 
within  5  percent  of  in  situ  values  are  obtained  (Ballard  and  McGee  199 1, 
McGee  and  Caulfield  1991). 

A  plot  of  the  impedance  function  versus  laboratory  measurements  of 
density  from  core  samples  taken  in  the  Mississippi  Sound  is  presented  in 
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Figure  1.  Hamilton's  data  (represented  by  the  solid  line  in  Figure  1),  al¬ 
though  obtained  along  the  shelf  and  slope  of  the  continental  terrace,  show 
remarkable  agreement. 

At  this  point  if  should  be  noted  that  none  of  the  above  data  were  pro¬ 
cessed  by  matched  filter  correlation  procedures.  Marked  improvements  in 
data  resolution  have  been  noted  since  signal  correlation  has  been  im¬ 
plemented.  This  subject  will  be  discussed  in  more  detail  later. 

Case  History 

During  summer  1992,  an  acoustic  impedance  survey  of  Galveston  Ship 
Channel  was  conducted.  Although  several  objectives  were  accomplished, 
the  phase  concerning  volumetric  determinations  of  "dirty"  sands  (density 
1.7  to  1.95  g/cm3)  underlying  the  proposed  channel  were  of  specific  inter¬ 
est  because  of  their  potential  use  for  beaches.  The  following  discussion 
presents  the  methodology  used  to  locate,  identify,  and  quantify  that  mate- 
riak  of  interest  to  the  U.S.  Army  Engineer  District,  Galveston. 

Data  Acquisition 

In  the  normal  Course  of  data  acquisition,  field  records  related  to  ampli¬ 
tude  of  recorded  signal,  time,  and  distance  provide  the  geophysicist/ 
engineer  with  quick-look  assessments  of  data  quality  and  subbottom  condi¬ 
tions.  Data  are  oftentimes  dually  recorded  by  analog  and  digital  systems. 
Analog  presentations  are  usually  in  shades  of  gray,  while  digital  data  arc 


displayed  in  color.  Figure  2  is  a  typical  3.5-kHz  pinger  amplitude  cross 
section  obtained  in  Galveston  Ship  Channel.  Note  that  the  top  of  the 
graph  is  not  the  water  surface,  but  an  assigned  water  column  delay.  This 
offset  allows  full  vertical  expansion  of  the  subbottom  display,  which  in 
this  case  extends  into  the  subbottom  more  than  40  ft.  Changes  in  strati¬ 
graphy  are  readily  apparent.  Records  of  this  type  are  used  as  quick-look 
guidance  in  boring  placement. 

Upon  determination  of  the  reflection  coefficients  and  impedance  func¬ 
tion  at  known  locations,  the  virtually  continuous  seismic  profiles  are 
processed.  The  single-channel,  digitally  recorded  data  are  read  into  the 
processing  software  developed  with  Caulfield  Engineering  (Caulfield 
1991b)  and  corrected  for  transmission  losses  due  to  spherical  spreading 
and  compensated  for  absorption  losses  in  each  layer.  Further  data  en¬ 
hancement  is  accomplished  by  use  of  recently  developed  (Caulfield  1992) 
matched  filter  signal  correlation  processing.  Use  of  the  new  correlation 
procedures  has  resulted  in  signal-to-noise  ratio  improvements  of  at  least 
10  db. 


Correlation  processing  also  allows  for  identification  of  the  frequency 
characteristics  of  source  signals  and  noise,  and  provides  techniques  to  opti¬ 
mize  source  and  array  selection,  thus  improving  frequency  penetration  re¬ 
covery.  This  has  resulted  in  vertical  resolution  improvement,  reduced 
noise  contamination,  and  improved  ability  to  target  special  materials  or 


Classical  multilayer  algorithms  are  used  to  compute  equivalent  reflec¬ 
tion  coefficients  and  impedances  along  the  profile.  This  in  turn  provides 
density  estimates  of  shallow  subbottom  layers  and  classifies  the  litho- 
stratigraphy  (Hamilton  1972,  Caulfield  and  Yim  1983).  The  results  are  cor¬ 
rected  for  tidal  fluctuations  and  correlated  with  survey  positioning  data. 
Processed  results  are  presented  in  the  form  of  annotated  amplitude  cross 
sections  or  two-  and  three-dimensional  (3-D)  views,  color-coded  according 
to  material  density. 

Density  Predictions 

By  incorporating  the  virtually  continuous  coverage  of  subbottom 
materials  with  digital  terrain  modeling  techniques,  rapid  and  accurate 
computations  can  be  made  of  volume  and  material  type  to  be  removed  by 
dredging.  Furthermore,  a  detailed  database  has  now  been  established  for 
project  monitoring  and  long-term  planning.  Computed  sediment  densities 
within  the  project  area  can  be  displayed  in  a  color-coded,  3-D  view  as 
shown  in  Figure  3,  if  desired  by  the  user.  In  this  example,  lighter  shad¬ 
ings  are  indicative  of  less  dense  material;  the  darkest  are  analogous  to  den¬ 
sities  1.7  g/cirri. 

Displays  of  this  type  provide  much-improved  data  interpretation  and 
visualization  for  the  end  product  user  as  compared  to  standard  two- 
dimensional  presentations  generated  exclusively  from  boring  information. 
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However,  caution  must  be  exercised  by  the  user  to  maintain  an  awareness 
of  real  data  versus  computer-generated  extrapolations. 

Volumetric  Calculations 


Before  computer-assisted  volume  estimates  can  be  calculated,  a  continu¬ 
ous  3-D  computer  model  of  the  subbottom  data  must  be  generated  for 
each  survey  line.  In  addition,  a  3-D  perspective  model  (Figure  3)  consist¬ 
ing  of  a  composite  of  data  from  ail  individual  survey  lines  may  be  created 
for  use  in  modeling  proposed  channel  cuts,  evaluating  slope  stability,  and 
so  forth.  At  the  project  planner's  discretion,  he  may  elect  to  view  an  area 
of  interest  from  various  angles  or  create  different  displays  by  stripping  or 
slicing  at  any  desired  coordinate. 


Using  Figure  3  for  our  example,  the  volume  of  any  material  to  be  re¬ 
moved  can  be  easily  calculated.  This  example  will  predict  the  configu¬ 
ration  of  silty  ("dirty")  sands  (>1.7  g/cm3  density;  2,864  to  3,052  x  102 
g/cm2/sec  acoustic  impedance)  underlying  a  selected  segment  of  the  proj¬ 
ect  study  area.  Calculating  the  volume  of  material  present  within  a  se¬ 
lected  area  of  the  perspective  model  is  accomplished  by  calculating  the 
volume  of  material  present  within  the  corresponding  area  of  each  profile 
line  model. 


In  our  example,  the  profile  line  model  shown  in  Figure  4  corresponds 
to  the  section  of  the  3-D  perspective  model  (Figure  3)  nearest  the  viewer, 
Note  that  the  axes  labels  in  Figure  3  show  northing  and  easting  coordi¬ 
nates,  while  the  axis  label  in  Figure  4  reflects  relative  feet  from  start-of- 
line.  This  conversion  is  necessary  to  correlate  the  profile  line  models  to 
amplitude  records  that  are  plotted  with  locations  relative  to  start-of-Iine. 


Before  calculating  volumes,  the  area  of  interest  must  be  sliced  out  of  the 
computer  model  and  the  material  density  range  to  be  displayed  selected. 
In  this  example,  the  volume  of  material  present  in  the  area  between  the 
bottom  to  “55*0=ft  depth,  and  between  3,000  and  7,000  ft  from  st  art-o  f-line  f 
having  a  density  of  >1.7  g/cm3,  will  be  selected.  The  section  of  the 
model  below  -55  ft  is  stripped  away,  and  the  sections  from  0  to  3,000  ft 
and  from  7,000  to  19,000  ft  would  be  sliced  from  the  display,  as  shown  in 
Fieurc  5.  1 


Finally,  the  density  of  material  to  be  displayed  would  be  set  to  >1.7 
g/cm3.  Calculations  based  on  these  parameters  are  then  computer- 
generated  and  displayed,  as  shown  in  Figure  6,  yielding  an  estimated  vol¬ 
ume  of  167,366  cu  yd  of  "dirty"  sand  (materials  >1.7  g/cm1)  within  the 
specified  boundaries.  This  step  is  repeated  for  each  profile  line  model 
within  the  area  of  interest,  summing  the  estimated  volumes  for  an  overall 
project  volume  total. 
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Conclusions 

In  its  present  state  of  development,  acoustic  impedance  processing  of 
seismic  reflection  data  provides  an  accurate,  continuous  description  of 
bottom  and  subbottom  marine  sediment  characteristics  in  a  rapid, 
cost-effective  manner.  Results  from  properly  calibrated  surveys  have  been 
used  to  provide  Corps  Districts  and  dredging  contractors  with 

•  Density  estimates  of  marine  sediments. 

•  Continuous  subbottom  information  for  planning  and  designing  dredg¬ 
ing  and  sampling  programs. 

•  Estimates  of  the  volume  and  type  of  material  to  be  removed  through 
dredging. 

•  A  detailed  and  continuous  geologic  database  for  aiding  long-term  plan¬ 
ning  of  future  work. 

Acoustic  impedance  information,  if  properly  implemented  in  the  project 
planning  stages,  provides  valuable  data  on  the  distribution  and  extent  of 
differing  marine  sediments,  aids  in  locating  optimal  placements  of  sam¬ 
pling  cores,  and  supplements  previously  obtained  soil  borings  by  provid¬ 
ing  continuous  profile  coverage  of  sediment  characteristics  between 
sample  locations. 

Epilogue 

i  .  1  '  '  .  1  :  •  ■  ■  .  •  •  .  V.  f 

During  preparation  of  this  technical  note,  plans  were  formulated  to  per¬ 
form  high-resolution  acoustic  impedance  surveys  prior  to  dredging  opera¬ 
tions  in  Boston  Harbor,  Baltimore  Harbor,  and  Houston  Ship  Channel. 
Additionally,  an  acoustic  survey  is  being  planned  to  locate  gravel  beds  in 
the  Atchafalaya  River. 

Possibly  for  the  first  time  in  Corps  of  Engineers  history,  the  Mobile  Dis¬ 
trict  advertised  site  conditions  for  dredging  at  Culfport  Ship  Channel 
using  density  predictions  determined  by  geophysical  acoustic  impedance 
surveying  in  addition  to  traditional  bore  hole  information.  True  site 
conditions,  resulting  from  dredging  operations  begun  in  May  1992,  have 
been  closely  monitored  and  compared  to  predictions,  Results,  as  reported 
by  Mobile  District  inspectors  and  dredging  contractor  Bean,  have  been 
excellent. 

References 

Ballard,  R.  F.,  Jr.,  and  McGee,  R.  G.  1991.  "Subbottom  Site  Characterization 

by  Acoustic  Impedance,"  Proceedings,  U.S.  Army  Corps  of  Engineers  Survey¬ 
ing  Conference,  15-17  July  1991,  Louisville,  KY. 


Technical  Note  DRP-2-07  (July  1993)  9 


Ballard,  R.  P.,  Jr.,  McGee,  R.  G.,  and  Whalm,  R.  W.  1992.  "A  High-Resolution 
Subbottom  Imaging  System,"  Proceedings ,  US/Japan  Marine  Facilities  Panel  of 
UJNR,  25  October-1 1  November  1992,  Washington,  DC. 

Caulfield,  D.  D.  1991a.  "Acoustic  Core  System  Manual,  Version  4.5x,"  Caul¬ 
field  Engineering  Group,  Oyama,  British  Columbia,  Canada. 

Caulfield,  D.  D.  1991b.  "Digital  Field  Shallow  Seismic  Acquisition  System 
Manual,  Version  1.2x,"  Caulfield  Engineering  Group,  Oyama,  British  Co¬ 
lumbia,  Canada. 

Caulfield,  D.  D.  1992.  "Final  Report — Improvement  in  Seismic  Resolution 
Using  Correlation  Techniques,"  Caulfield  Engineering  Group,  Oyama, 
British  Columbia,  Canada. 

Caulfield,  D.  D.,  and  Yim,  Y.  C.  1983.  "Prediction  of  Shallow  Subbottom  Sedi¬ 
ment  Acoustic  Impedance  While  Estimating  Absorption  and  Other 
Losses,"  Journal  of  the  Canadian  Society  of  Exploration  Geophysicists,  Vol  19, 

No.  1,  pp  44-50. 

Hamilton,  E.  L.  1970.  "Sound  Velocity  and  Related  Properties  of  Marine  Sedi¬ 
ments,  North  Pacific,"  Journal  of  Geophysical  Research ,  Vol  75,  No.  23,  pp 
4423-4446. 

Hamilton,  E.  L.  1972.  "Compressional-Wave  Attenuation  in  Marine  Sedi¬ 
ments,"  Geophysics,  Vol  37,  No.  4,  pp  620-646. 

Hamilton,  E.  L.,  and  Bachman,  R.  T.  1 982.  "Sound  Velocity  and  Related  Prop¬ 
erties  of  Marine  Sediments,"  Journal  of  Acoustics  Society  of  America,  Vol  72, 
No.  6,  pp  1891-1904. 

McGee,  R.  G.,  and  Caulfield,  D.  D.  1991.  "Practical  Implementation  of  the 
Acoustic  Core  Density  Estimating  System,"  1991  National  Convention, 
Canadian  Society  of  Exploration  Geophysicists,  Calgary,  Alberta,  Canada. 


10 


Technical  Note  DRP-2-07  (July  1993) 


Figure  2.  Pinger  (3.5-kHz)  amplitude  cross  section,  Galveston  Ship 

Channel  ' 


Figure  3.  Three-dimensional  perspective  view  of  inner  bar  and 

anchorage  basin,  Galveston  Harbor.  Axes  are  displayed  in  northing 
and  easting  coordinates 


Figure  4,  Full  profile  mode! 


Figure  6.  Profile  model  showing  only  material  1.7  g/cm3  density  within  area  of 

interest 
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